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Abstract 

Accounting for double-logarithms of x and running QCD coupling leads to 
expressions for both the non-singlet and singlet components of gi. These ex- 
pressions manifest the Regge asymptotics when x and differ considerably 
from the DGLAP expressions at small values of x. 



1 Introduction 



As is well known, deep inelastic scattering (DIS) is one of the basic processes 
for probing the structure of hadrons. From the theoretical point of view, the 
inclusive cross section of DIS is a convolution of the leptonic and hadronic 
tensors, with the information about the structure of the hadrons participating 
into DIS coming from the hadronic tensor . the forward Compton amplitude, 
when a deeply off-shell photon with virtuality scatters off an on-shcU hadron 
with momentum p. The spin-dependent part, of the hadronic tensor is 

parametrized in terms of two structure functions, 171 and 52, as 



^ui. =^^tLu\p *p.gi+(jp Pp)92 ~i^^iu\p ^pgi+^o [91+92] 

pq I- pq i pq I ^ ' 

(1) 

so that gi is related to the longitudinal hadron spin-flip scattering, whereas the 
sum gi + (72 is relative to the transverse spin- flips. In Eq. J^l, ™ stands for 
the hadron mass, Sp and are the longitudinal and transverse (with respect 
to the plane formed by p and q) components of the hadron spin Sp. Both gi 
and 52 depend on x = —q^/2pq, < a; < 1 and = —q"^ > 0. Obviously, 
small X corresponds to s = [p + q)^ ~ 2pq 3> Q^. In this case, Sp ~ Pp/m 
and therefore the part of W^J^ related to g\ does not depend on m. When 
^ m^, one can assume the factorization and regard as a convolution 

of two objects (see Fig. 1). The first one is the probability <&(<!> = $g in 
Fig. la and $ = <I>g in Fig. lb) to find a polarized parton (a quark or a 
gluon) within the hadron. The second one is the partonic tensor defined 

and parametrized similarly to W^^"\ Whereas the partonic tensor W^"\ i.e. 
the partonic structure functions g\ and (72, can be studied within perturbative 
QCD, ^q.g are essentially non-perturbative objects. The lack of knowledge of 
$ is usually compensated by introducing initial parton distributions which can 
be found from phenomenological considerations. In doing so, the non-singlet 
component g^^ of g\ is usually expressed as a convolution of a piece which we 
denote g^^ from purely evolution, and the initial polarized quark density Ag: 

g^ = ^ ® Aq . (2) 

Similarly, the singlet structure function gf is expressed in terms of the evolution 
pieces gq and gg and the densities of the polarized quarks and gluons, Aq and 
Ag: 

5f ^gq(^Aq + gg ® Ag . (3) 

The subscripts q, g in Eqs. (|2I3|I refer to the kind of the initial partons, as is 
shown at Fig. 1. We remind that there is no rigorous procedure for calculating 




Figure 1: Representation of the hadronic gi as the convolution of the partonic 
gi and the initial parton densities. 



Ag and lS.g. They have to be found from phenomenological considerations. 
On the contrary, there are regular perturbative methods for calculating the 
evolution parts g^^ and gq^g of gi. When is much greater than the starting 
point of the Q^-evolution and at the same time a; ^ 1, it is convenient to 
rewrite Eqs. (|2I3|I in the form of the Mellin integral: 
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(4) 



with CNs{^)iCq^g{Lij) being the coefficient functions and 7Ars(w), 7s(w) the 
non-singiet and singlet anomalous dimensions respectively. The anomalous 
dimensions control the Q^-evolution and the coefficient functions govern the 
s-evolution which, at fixed , is equivalent to the x-evolution. 



The best known instrument to calculate the DIS structure functions is 
the DGLAP approach. In the DGLAP framework, both the coefficient func- 
tions and the anomalous dimensions are perturbatively known and represented 
by their one-loop, (or leading order (LO)) ^ and two-loop (next-to-leading 
order (NLO))contributions. The relevant Ref.s concerning the NLO calcula- 
tions can be found in the review El . The remaining ingredients to the rhs 
of Eq. igj), Ag and Ag can be taken, for example, from Ref. 0. DGLAP 

provides a quite good description of the experimental data"^. The extrap- 
olation of DGLAP into the small-x region predicts an asymptotic behavior 
^ exp(-\/Cln(l/a;) In In Q^) for all DIS structure functions (with different fac- 
tors C). However, from a theoretical point of view, such an extrapolation at 
the small- a; is rather doubtful. In particular, it neglects in a systematical way 
contributions ot the type ^ (as ln^(l/a;))'^ which are small when x ^ 1 but be- 
come relevant when a; <C 1. The total resummation of these double- logarithmic 

(DL) contributions was made in Refs. El and Ref. ^ for the non-singlet (gf ■^) 
and singlet gi respectively, and it leads to the Regge (power-like) asymptotics 

9i{9i^) ^ (l/a;)'^"''((l/a:;)'^"^), with A^^,A^^ being the intercepts calcu- 
lated ib the double- logarithmic approximation (DLA). The weak point of this 

resummation in Refs. El ^ El ig the assumption that is kept fixed (at some 
unknown scale) . It leads therefore to the value of the intercepts A^^ , A^^ 
explicitly depending on this unknown coupling, while is well-known to be 

running. The results of Refs. El^El had led many authors (see e.g. El) to sug- 
gest that the DGLAP parametrization = ctsiQ'^) has to be used. However, 
according to results of Ref. ^ such a parametrization is indeed correct for 
I ~ 1 only and cannot be used for cc <C 1. The appropriate dependence of as 
suggested in Ref. I3> , has been successfully used to calculate both and gi 
singlet at small x in Refs. El . in the present talk we review these results. 

Instead of a direct study of gi like it is done in DGLAP, it is more con- 
venient to consider the forward Compton amplitude A for the photon-parton 
scattering related to gi, with Aq = 1, Ag = 0, as follows: 

gi{x,Q^)^-^sA{s,Q') . (5) 

IT 

As already stated above, we cannot use DGLAP for studying gi or A at 
small x because it does not not include the total resummation of the double- and 
single-logarithms of x for the anomalous dimensions and coefficient functions, 
and also the as-parametrization used is valid for x not far from 1. 

Then, in order to account for the double-logs of both x and Q^, we need 
to construct a kind of two-dimensional evolution equations that would combine 
both the X- and Q^- evolutions. 



Such equations should sum up the contributions of the Feynman graphs 
involved to all orders in as- Some of those graphs have either ultraviolet or 
infrared (IR) divergences. The ultraviolet divergences are regulated by the 
usual renormalization procedure. In order to regulate the IR ones, we introduce 
an IR cut-off in the transverse momentum space for the momenta ki of all 
virtual quarks and gluons: 

M < h± (6) 

where ki± stands for the transverse (with respect to the plane formed by the 
external momenta p and q) component of ki . This technique of regulating the 

IR divergences was suggested by Lipatov and used first in Ref. ^ for quark- 
quark scattering. Using this cut-off /i, A acquires a dependence on /i. Then, one 
can evolve A with respect to /i, constructing the appropriate Infrared Evolution 
Equations (IREE). As A ^ A{s / fi'^ , / ^'^) , 

- fi^dA/dfi^ = dA/dp + dA/dy (7) 

where p = ln(s//i^) and y — \n{Q^/p^). Eq. j?)) represents the Ihs of the 
IREE for A. In order to write the rhs of the IREE, it is convenient to use the 
Sommerfeld- Watson transform 

/too 1 
— is/p^racj)F{cj,Q^) (8) 
-too 27ri 

where ^(w) is the negative signature factor, ^(w) = [1 — e^*'^"]/2 « nTLu/2. It 
must be noted that the transform inverse to Eq. JHJ involves the imaginary 
parts of A: 

F(w,Q2) = _/ dpe-P'^QAis^Q^) . (9) 

Notice that, contrary to the amplitude A, the structure function gi does not 
have any signature and therefore ^(w) = 1 when the transform (|SJ is applied 
directly to gi. 

2 infrared evolution equations for gi 

When the factorization depicted in Fig. 1 is assumed, the calculation of gi (for 
semplicity we will omit the superscript "s" iorgi singlet, though we use the 
notation gf^ for the non-singlet gi ) is reduced to calculating the Feynman 
graphs contributing the partonic tensor Wj^^ depicted as the upper blobs in 
Fig. 1. Both cases, (a), when the virtual photon scatters off the nearly on-shell 
polarized quark, and (b), when the quark is replaced by the polarized gluon. 



should be taken into account. Therefore, in contrast to Eq. 0, we need to 
introduce two Compton amphtudes: Aq and Ag corresponding to the upper 
blob in Fig. la and Fig. lb respectively. The subscripts "q" and "g" refer to 
the initial partons. Therefore, 



(10) 



1 



where 

gq = -5,A,(s, Q^), gg = -5,Ag(s, Q') (11) 

TT TT 

Let us now construct the IREE for the amplitudes Aq^g related to gi. To 
this aim, let us consider a virtual parton with minimal /cj_. We call such a 
parton the softest one. If it is a gluon, its DL contribution can be factorized, 
i.e. its DL contribution comes from the graphs where its propagator is attached 
to the external lines. As the gluon propagator cannot be attached to photons, 
this case is absent in IREE for Aq g. The second option is when the softest 
partons are a t-channel quark-antiquark or gluon pair. It leads us to the IREE 
depicted in Fig. 2. Applying the operator —^^d/dfi^ to it, combining the result 
with Eq. (O and using JSJ, we arrive at the following system of equations: 



d 



-^[Fqq(u;)Fq{LO,v)+Fqg{Lo)Fg{u;,y)\ 



Fgg{cj)Fg{u;,y)] . (12) 



The amplitudes Fq,Fg are related to Aq,Ag through the transform (|SJ|. The 
Mellin amplitudes Fik, with i,k = q,g, describe the parton-parton forward 
scattering. They contain DL contributions to all orders in as- We can intro- 
duce the new anomalous dimensions Hn^ — (l/87r^)i^ifc. The subscripts "q,g" 



correspond to the DGLAP-notation. 
Eq. (^J leads to 



Solving this system of eqs. and using 
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The unknown factors C±{lu) have to be specified and will be discussed 
later. All other factors in Eq. (|13l) can be expressed in therms of Hn,: 



X — Hgg Hqq , 



^ - (Hgg - Hqqf + ^HqgHgq , 



(14) 
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Figure 2: Infrared evolution equations for the amplitudes Aq,Ag. 



The anomalous dimension matrix Hih was calculated in Ref. 
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We have used here the notations Cp = 4/3, iV = 3 and Uf = A. The quantities 
A(uj) and D(uj) account for the running of as. They are given by the foUowing 
expressions: 



A{u;) 
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dpe 



rf + TT^ 7q [p + 77)2 + TT^ 



P + V 



{p + riY + ir'^ p + T] 



(21) 
(22) 



with T] = \n{p? / Aq^jj) and b — (33 — 2rij)/127r. A' is defined as A with the 
TT^ term dropped out. 

FinaUy we have to specify the coefficients functions C± appearing in 
Eq. (El. When = p^ , 



9g = A5(xo) 



(23) 



where Ag(a::o) a-nd Ag{xQ) are the input distributions of the polarized partons 
at xq = /i^/s. They do not depend on Q^. Eq. (|23() aUows us to express C±{uj) 
in terms of Aq{uj) and Ag{Lu), which are related to Ag(a;o) and Ag{xo) through 
the ordinary Mellin transform. Indeed, 



c+ + C- = Aq, a 
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2i/„ 
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(24) 
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with both Aq and Ag depending on uj. This leads to the following expressions 
for gq and gg-. 

ggix,Q^)= |^(l/a;)-[(^(-)Ag+i?A5)e"+^+(A(+)Ag-BA5)e^-^] , 

(25) 



with 



/ —(l/xr\(EAq+A^+^Ag)e^+y+(-EAq+A^-^Ag 
V2 2\/i?/ Vi? VR 



(26) 
(27) 



Eqs. H25I (|26|l express 171 in terms of the parton distributions Aq{uj) and 
Ag{uj), which are related to the distributions Aq{xQ) and Ag{xo) at very low x: 
Xq « n"^ / s <C 1. Therefore, they hardly can be found from experimental data. 
It is much more useful to express gq , gg in terms of the initial parton densities Sq 
and Sg defined at a; '--^ 1. We can do it, using the evolution of Aq{xQ), Ag{xo) 
with respect to s. Indeed, the s-evolution of Sq,Sq from s ~ /i^ to s ^ /i^ at 
fixed [Q^ = M^) is equivalent to their x-evolution from a; ~ I to x <C 1. 
In the cj-space, the system of IREE for the parton distributions looks quite 
similar to Eqs. (|12|l . However, the eqs for Aq, Ag are now algebraic because 
they do not depend on Q^: 



Aq{Lu) = {<el> /2)5q{uj) + [l/u) [Hgq{uj)Aq{uj) + Hqg{uj)Ag{oj)] , 
Ag{Lo) = i<el> /2)6g{uo) + {l/uj) [Hgq{uj)Aq{uo) + Hgg{u)Ag{uj)] (28) 

where < > is the sum of the quark electric charges (< >= 10/9 for 

nf = 4) and Sg = —{A'{u!)/TTUj'^)Sg is the starting point of the evolution of the 
gluon density 6g. It corresponds to Fig. lb where the upper blob is substituted 
by the quark box. Solving Eqs. H28|) . we obtain: 

^ (< eg > /2) HcJ -~ Hgg)Sq + cjHqgSg] 

[u;^-U;{Hgq + Hgg) + {HggHgg~HggHgq)] ' 

^ (< eg > /2) [ujHggSq + cjju; - Hgg)Sg] ^^^^ 

- U;{Hqq + Hgg) + [HqgHgg - HqgHgq)] 

Then Eqs. H25I2()I29I3U|I express gi in terms of the initial parton densities 
5q,5g. 

When we put Hgg = Hgq — Hgg = and do not sum over Cg, we arrive 
at the expression for the non-singlet structure function g^^: Obviously, in this 
case v4(+) = B = E = ri_ = 0, A'-~^ = 1, — Hqq. However, the nonsinglet 
anomalous dimension Hqq should be calculated in the limit hgg = hqg — bgq = 0. 
We denote such Hqg = H^^ . The explicit expression for it is: 



if^^ = (1/2) u; - Joj^ - Abqg . (31) 



Therefore, we arrive at 
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,^5 
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Figure 3: Dependence of the leading singularity uiq on ij. 



3 Small- a; asymptotics for gi 

Before making use of Eqs. (|25|) . H26l and (|32|l for calculating gi at small but finite 
values of x, let us discuss its asymptotics. When a; ^ and » /z^, one can 
neglect contributions with f2_ in Eqs. (|13|l . As is well known, gi ^ at 
a; ^ 0, with ujq being the position of the leading singularity of the integrand 
of gi . According to Eqs. (|15|l . the leading singularity, cu^^ for g^^ is the 
rightmost root of the equation 

a;2 _ ^i,^^ ^ (33) 
while the leading singularity, ujq for gi is the rightmost root of 

LO^ - 4(6,, + bgg)LU^ + I6{bgqbgg - bggbgg) = . (34) 

In our approach, all factors bik depend on i] — hi(^'^ / Aqcd), so the roots 
of Eqs. H33I34|I also depend on r/. This dependence is plotted in Fig. 3 for ui^^ 
and in Fig. 4 for luq. Both the curve in Fig. 3 and the curve 1 in Fig. 4 have 
a maximum. We denote this maximum as the intercept. Therefore, 

5f^~e^5g(l/x)^-(gVA^')^"^/^ (35) 

and we find for the intercepts 

Ans ~ 0.4, As ~ 0.86 (36) 
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Figure 4: Dependence of the leading singularity luq on rj. Curve 2 corresponds 
to the case where gluon contributions only are taken into account; curve 1 is 
the result of accounting for both gluon and quark contributions. 



and Zi = —1-2, Z2 ~ —0.08. This implies that Qi is positive when x 
whereas gf can be either positive or negative, depending on the relation 
between 5q and 5g. In particular, gi is positive when 

lb5q + 6g<Q. (37) 

otherwise it is negative. In other words, the sign of 51 at small x can be positive 
if the initial gluon density is negative and large. 



4 g^ at finite values of x 

Let us estimate the impact of the total resummation of DL and SL contributions 
on g^^. According to Eq. (|32|l . the value of gi{x,Q^)^^ depends both on 
the perturbative terms and on the inputs Aq. The latter can be obtained by 
fitting the experimental data and it is known (see e.g. ^) that widely different 
formulae for Aq can be used. In order to avoid discussing the fitting procedure 
and as in this paper we present our results only the perturbative part of g^^ , 
we can assume 

Aq = S{1 - x) . (38) 
Then let us calculate the ratio 



R = , (39) 

where gi'^^ is the LO DGLAP non-singlet gi. The results of a numerical 
calculations for R aX — 20GeV^, ji = 1.5GeV are shown in Table 1. 



Table 1: Ratio R of structure functions g^^ calculated in DLA and DGLAP 
LO. 



X 


R 


0.1 


1.60 


0.01 


2.59 


0.001 


4.33 


0.0001 


7.46 



It shows that the hiipact of the total resummation of DL contributions 
is negligible for a; > 0.1 but it grows fast with decreasing of x, achieving the 
values i? = 2.6 at x = 10"^ and i? = 7.5 at x = 10"''. 

5 Conclusion 

The total resummation of the most singular (~ Q!"/(jj^"+^) terms in the expres- 
sions for the anomalous dimensions and the coefhcient functions leads to the 
expressions of Eqs. (|10I25I26I32|I for the singlet and the non-singlet structure 
functions gi. It guarantees the Regge (power-like) behavior (|35|l of gi, g^^ 
when X —> 0, with the intercepts given by Eq. H36(l . The intercepts A^rs, A5 
are obtained with the running QCD coupling effects taken into account. The 
value of the non-singlet intercept A^s ~ 0-4 is now confirmed by several in- 
dependent analysis of experimental data and our result A5 w 0.86 is in a 
good agreement with the estimate of Ref. : Ag = 0.88 ±0.14 obtained from 
analysis of the HERMES data. Eq. states that gf^ is positive both at 
X ^ 1 and at a; ^ 1. The situation concerning the singlet gi is more involved: 
being positive at a; ~ 1, the singlet gi can remain positive at x <^ 1 only if the 
initial parton densities obey Eq. (|S7|l . otherwise it becomes negative. The ratio 
of our results versus the DGLAP ones for non-singlet gi is given in Table 1. 
It shows explicitly that the impact of high-order DL contributions is small at 
X > 0.1 but it grows fast when x is approaching 10^"^ — 10~^. 
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